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We analyze the recent experiment [W. Neubeck et al., Phys. Rev. B 60 (1999) R9912] for the 
resonant x-ray magnetic scattering (RXMS) around the K edge of Co in the antiferromagnet 
CoO. We propose a mechanism of the RXMS to make the 4p states couple to the magnetic 
order: the intraatomic exchange interaction between the 4p and the 3d states and the p-d 
mixing to the 3d states of neighboring Co atoms. These couplings induce the orbital moment 
in the 4p states and make the scattering tensor antisymmetric. Using a cluster model, we 
demonstrate that this modification gives rise to a large RXMS intensity in the dipole process, 
in good agreement with the experiment. We also find that the pre-edge peak is generated by 
the transition to the 3d states in the quadrupole process, with negligible contribution of the 
dipole process. We also discuss the azimuthal angle dependence of the intensity. 
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§1. Introduction 

Resonant x-ray scattering has recently attracted much interest as a useful tool to investigate the 
magnetic order. The resonant enhancement for magnetic Bragg reflections has been observed in 
transition metalsS^ and their compounds.i'BB For the resonant enhancement to be observable on 
the magnetic Bragg spots, orbitals have to couple with the magnetic order. It is usually achieved 
through the spin-orbit interaction. In this paper, we study the mechanism of the resonant x-ray 
magnetic scattering (RXMS) by focusing our attention on a typical material, CoO, where orbitals 
are strongly coupled to the magnetic order, and analyze the recent RXMS experiment by Neubeck 
et a$ 

In CoO, Co atoms form an fee lattice with a small tetragonal distortion of c/a < 1 below the 
Neel temperature Tjy = 292 K. It is a type II antiferromagnet (AF) characterized by a wavevector 
Q, which is directed to one of four body-diagonals in the fee lattice, thus forming K domains.!'!) 
Although a non-collinear magnetic order has been proposed so far,@'& it is quite reasonable to 
assume a collinear magnetic order in the analysis of the RXMS intensity, since the RXMS signal 
is separately observed for different domains in the experiment .B The spin vector lies in the plane 



l 



perpendicular to the AF modulation direction; three directions are possible, forming S domains.^) 
For an example, the spins are directed to (—1,-1,2), (—1,2,-1), or (2,-1,-1), which we call 
domains Si, S2, and S3, in the K domain with Q = (5,5,5)- Figure 1 depicts schematically this 
domain. In total, we have 12 domains. The orbital moment is known to be as large as ~ 1/^_b-& 
This material was expensively studied for decades ago by Kanamori,!) who pointed out that the 
large orbital moment causes a large magnetostriction. 
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Fig. 1. Sketch of a CoO crystal (only Co atoms are shown). Wave vector Q = (\,\,\) characterizes an antifer- 
romagnetic modulation direction. Three arrows in a plane perpendicular to Q indicate possible spin directions. 
Spins of solid, gray, and open circles align alternately. 



For photons with energy near the K edge, the Ap states become occupied in the intermediate 
states of the dipole process. The 4p states, different from the 3d states, are highly extended in 
space. To have a feeling of the character of the wavefunctions, we show in Fig. 2 the radial wave 
functions calculated within the Hartree-Fock approximation© for a free Co 2+ atom. The concept 
of localized orbitals loses its clear meaning outside the ion radius in solids. Since the 4p states have 
large weights outside the ion radius, treating them as local orbitals in a cluster model may limit the 
usefulness of the model to a qualitative level. Nevertheless, we use a cluster model to make clear 
the mechanism of the RXMS in CoO. We hope this will serve as a first step to more quantitative 
study based on band structure calculations. 
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A closely related phenomenon where the Ap states are involved is the resonant x-ray scattering 
for the orbital order in LaMn0 3 fl@) This case has been thoroughly studied with treating the Ap 
states of Mn as a band-HHS© It has been shown that the Ap states of Mn is mainly modified by the 
Jahn- Teller distortion via the oxygen potential on the neighboring sites, and that this modification 
gives rise to non-vanishing diagonal elements in the scattering tensor, leading to a large main peak. 
On the other hand, the mechanism of the RXMS for CoO is much complicated. The Ap states are 
required to couple to the magnetic order. This is achieved through the coupling to the 3d states 
which constitute the magnetic order. Two processes are likely to be effective to make the Ap states 
couple to the 3d states: the intraatomic exchange interaction and the p-d mixing to the 3d states of 
neighboring Co atoms. We demonstrate that these couplings give rise to a modification of the Ap 
states and make the scattering tensor antisymmetric. The RXMS intensity thus generated seems 
much smaller than the intensity generated by the Jahn- Teller effect in LaMnOs, but is sufficiently 
larger than the pre-edge intensity coming from the quadrupole process. Note that the 3d states 
having orbital moment is crucial to generate the RXMS intensity. The amplitude of exciting an 
electron with up spin (up spin process) is generally different from that of the exciting an electron 
with down spin (down spin process). However, the sum of the amplitudes for the up and down spin 
processes become the same at each site if the system has no orbital moment; the RXMS intensity 
is not generated by such spin scattering processes alone, since the total amplitude is given by the 
difference between those at the A and B sublattices. 
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Fig. 2. Radial wave functions (multiplied by the radius) in the 3d 7 configuration of a free Co 2+ atom. The radius 
r is measured in units of the Bohr radius ao. The arrow a indicates the ion radius, which is defined by the half of 
the distance between the Co and O atoms in the CoO crystal. The arrow b indicates the distance to a neighboring 
O atom. 



In addition to the dipole process, we also consider the quadrupole process, which gives rise to a 
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pre-edge peak. Although the dipole process can also contribute to the pre-edge peak through the 
p-symmetric states mixing to the 3d states of neighboring Co atoms, this contribution is found to be 
much smaller than that of the quadrupole process. This result is consistent with the experiment 
Note that, for the pre-edge peak in LaMnOs, the contribution from the dipole process is predicted 
to be much larger than that of the quadrupole process© 

This paper is organized as follows. In § 2, we formulate the RXMS intensity. In § 3, we present 
a cluster model and discuss the relevant Ap and 3d states. In § 4, we show the calculated results 
for the RXMS intensity in comparison with the experiment. Section 5 is devoted to concluding 
remarks. In the Appendix, the geometrical factor is explicitly given for the dipole and quadrupole 
processes. 

§2. X-Ray Magnetic Scattering 

The magnetic scattering amplitude has been reviewed by several authors. We summarize 

here the expression of the amplitude in a form suitable to calculating the azimuthal angle depen- 
dence. The scattering geometry is shown in Fig. 3. Photon with frequency uj, momentum kj and 
polarization fj, (= a or ir) is scattered into the state with momentum kj and polarization // (= a' 
or it'). For the system with two magnetic sublattices A and B in the AF order, we write down the 
cross section for elastic scattering at magnetic Bragg peaks as 



da 
~dVL' 



oc 



T^ f (G,u)-T^r(G,u;)\ (2.1) 



with 



Z^(G,oO = - % QL"(G) • A" + S"(G) • B 

+ ^(G, W ) + L^(G >W ). (2.2) 

where G = kj — kj is the scattering vector. The first term in eq. (2.2) represents a non-resonant 
term at sublattice rj (= A or B), with m and c being the electron mass and the velocity of photon, 
respectively.0 L r? (G) and S 7? (G) stand for the Fourier transform of the orbital and spin angular 
momentum density,0 and A" and B are given by 

A" = A' — (A' • G)G, A' = -4 sin 2 9{e x e), (2.3) 
B = e x e + (k/ x e')(k/ • e) - (kj x c)(kj • e') - (k/ x e) x (kj x e), (2.4) 

where e and e are the initial and scattered polarizations, and kj, kj , and G are normalized vectors 
of kj, kj, and G. 

The second and third terms in eq. (2.2) represent the resonant contribution at the K edge, 
which originates from a second-order process that a photon is virtually absorbed by exciting the 
Is electron to unoccupied states and then emitted by recombining the excited electron with the Is 
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Fig. 3. Geometry of x ray scattering. Photon with polarization a or tt is scattered into the state of polarization a' or 
7r' at Bragg angle 8. The sample crystal is rotated by azimuthal angle ip around the scattering vector G = k/ — ki. 



core hole. The amplitude by the dipole process is expressed as 



(2.5) 



where the geometrical factors P^, P'^' are explicitly given in the Appendix, and the scattering 
tensor M v (ui) is given by 



**v f \ muJ jg(9\x a \j){j\x a '\g) 



-^Tioj- (Ej -E g ) + iT/2' 



(2.6) 



Here x a 's stand for the dipole operators, x\ = x, %2 = y, and x% = z in the coordinate frame 
fixed to the crystal axes. The ground state \g) has an energy E g , and the intermediate state \j), 
consisting of an excited electron on the Ap states and a hole on the Is state, has an energy Ej\ 
ujjg = (Ej — E g )/h. Quantity T is introduced to describe the life-time broadening due to the core 
hole. The amplitude by the quadrupole process is expressed as 

77' 

where the geometrical factors, Q^, Q'^\ are explicitly given in the Appendix, and the scattering 
tensor N v (u) is given by 



k 2 



zy\g) 



Huj - (E t - E g )+iT/2' 



(2.8) 



Here k is the wavenumber of the incident (and scattered) photon, and z^s stand for the quadrupole 
operators, z\ = (\/~3/2)(x 2 — y 2 ), Z2 = (1/2) (3z 2 — r 2 ), z% = V3yz, z^ = y/3zx, and z^ = V3xy in 
the coordinate frame fixed to the crystal axes. The intermediate state \£), consisting of an excited 
electron on the 3d states and a hole on the Is states, has an energy Ef, uj^ g = (Ef — E g )/h. 
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§3. 3d and 4p States at Co Sites 

Each Co atom is in the 3d 7 configuration under the cubic crystal field. We neglect small effects 
of the distortion such that c/a < 1. The ground state 4 F of a free cobaltous ion is split into the 
representations 1?4 and T§. Three degenerate states in the T4 representation have lower energy than 
those in the T§ representation. The spin of each ion is subjected to a strong exchange field and is 
oriented in its direction, forming an AF long range order. As mentioned before, we assume that 
the AF modulation direction Q is one of the four body-diagonals in the fee lattice. 

As regards the orbitals, we introduce the pseudo-angular-momentum operator I of magnitude 1 
in the T4 representation, such that the orbital-angular-momentum operator L is given by —(3/2)£. 
Thereby, in the first approximation, the wave function of the state specified by l z = correspond 
to the configuration that </?q = d xy is singly occupied (in other wards, d zx and d yz are doubly 
occupied); the state l z = ±1 is such that ip±\ == =f-^(d zx ± id yz ) is singly occupied. Connected 
with the spin through the spin-orbit coupling, I becomes antiparallel to the spin direction. Thus 
the ground state gives the spin moment ~ 3\ib and the orbital moment ~ 1.5//^. As discussed by 
Kanamori,& quantum mechanical effects modify this classical ground state, making the magnetic 
moment close to the experimental value. Since we are aiming at analyzing qualitatively the RXMS 
intensity, we will not consider the quantum correction. 

The unoccupied 3d states become occupied by an electron excited from the Is state in the 
intermediate states of the quadrupole process. The scattering amplitudes of filling the unoccupied 
3d states are usually cancelled out between those at the A and B sublattices. However, it is not the 
case for the unoccupied state in the T4 representation, since this orbital is strongly coupled to the 
magnetic order through the spin-orbit coupling. We call this unoccupied state in the intermediate 
state as the 3d* state with energy e^. This 3d* state, which is represented as in the coordinate 
frame that the z axis is along the spin direction (the upper (lower) sign corresponds to the A (B) 
sublattice), is expressed in the coordinate frame fixed to the crystal axes as 



where ,(a, P,-y) is the rotation matrix of the first rankj^ with the Euler angles 



(3.1) 



arcsm 



1 



a = — 37r/4, 

7 = for domain S\, a = tt/2 + arcsin 1/y/E, [3 = it — arcsin ^5/6, 7 = for domain 
S2, and a = — arcsin l/^/b, j3 = it — arcsin \/5/6, 7 = for domain S3. Substituting these states 
for \£) in eq. (|2.8D, we obtain the scattering tensor, which has an antisymmetric form, 



N (u>) — N B (ui) 
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with d, e, and / being complex numbers. 

The unoccupied Ap states at Co sites become occupied by an electron excited from the Is state 
in the intermediate states of the dipole process. The scattering amplitudes of filling the Ap states 
are cancelled out between those at the A and B sublattices, unless the Ap states are coupled to the 
3d states which constitute the magnetic order. One possible coupling is the intraatomic exchange 
interaction between the Ap and the 3d states. Since the occupied 3d states have finite orbital 
moment, this coupling makes the Ap levels split; at the A sublattice, the states p[ = ~-^(p' x ~\~Wy): 
Po = Pz an d P-i = ~^(p'x ~ iPy) with down spin have energy e\, 0, 62] at the B sublattice, p[, p' , 
p'_ 1 with up spin have energy €2, 0, e\. Here p'±± and p' refer to the coordinate frame that the z axis 
is along the spin direction. The ei and €2 are roughly evaluated as e\ = —0.056 eV, €2 = —0.020 
eV.© Another coupling is the mixing to the 3d states of neighboring Co atoms. What is effective 
to make the amplitude at the A sublattice different from at the B sublattice is the mixing to the 
3d* state. 

To estimate the effects of these couplings, we use a cluster model that the Ap states at the central 
site are mixing to the 3d* state of 12 neighboring Co atoms (Fig. 1). Thus the Hamiltonian matrix 
has 39 x 39 dimensions. The effect of the mixing to the 2p states of neighboring oxygens is implicitly 
included into the model by adjusting the energy of the Ap states. According to the experiment,!^ 
the main peak is at huj = 7.724 keV and the pre-edge peak at Tluj = 7.707 keV. Thereby we put 
£p — £d = 17 eV for the difference between the energy of the Ap states, e p , and that of the 3d* states, 
ed, both at the central site. Note that both e p and are including the interaction energy with the 
Is hole. The energy of the 3d* state at the neighboring Co atoms is several eV's higher than that 
at the central site where the core hole exists, since the interaction with the Is hole is not working 
there. We assume the difference to be Id — e<i = 4 eV. Then we evaluate the off-diagonal part of the 
Hamiltonian matrix by using the Slater-Koster parameters, (pda) and (pdyr)^ with the help of 
eq. (|3.1| ). We assume (pda) = —0.5 eV together with an empirical relation (pdir) = — (1/2) (pdcr)© 
This size of the (pda) value seems reasonable, since the Ap states are highly extended in space. 

The Hamiltonian matrix can be numerically diagonalized. We find that three degenerate Ap 
states are split into the states with finite orbital moments. Substituting these eigenstates for the 



intermediate states in eq. (2.6), we calculate the scattering amplitude. As a reflection of the induced 
orbital moment for the Ap states, the scattering tensor has an antisymmetric form, 



(3.3) 
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with a, b, and c being complex numbers. 
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§4. Calculated Results in Comparison with Experiment 

In the preceding section, we have studied the 3d and 4p states in the intermediate states of 
the RXMS. For calculating the scattering intensity, we have to estimate further several quantities. 
First, we evaluate the matrix elements for the dipole and quadrupole processes within the HF 
approximation© They are given by 



where Ri s (r), i?4 p (r), and i?3^(r) are radial wave functions for the Is, 4p, and 3d states in the 
3d 7 configuration of a free cobaltous ion (see Fig. 2). Since the Is state is well localized inside the 
ion radius, these values are to be good estimates in solids, with multiplying a factor 0.8. With 
these values and fc~3.9x 10 8 cnr 1 around the K edge, we have /c 2 |(3d|r 2 |ls)| 2 /(20|(4p|r|ls)| 2 ) ~ 
1.5 x 10 -2 , which is a measure of the ratio of the amplitude of the quadrupole process to that of the 
dipole process. Second, we assume that |S(G)| ~ 1.0, |L(G)| ~ 0.8, in the first term of eq. (2.3), 
for the scattering vectors G = (— i, — |) and (| , \ , |); these values are smaller than those for 
G = 0, that is, the local orbital- and spin-angular momenta. This is consistent with the general 
tendency that |S(G)| and |L(G)| are decreasing with increasing values of G@> The precise values 
are not necessary for our qualitative analysis. Third, we use the experimental value 1.7 eV from 
the fluorescence measurement for the life-time broadening width r.@) 

Experimentally, Neubeck et al. have separated the contributions to the scattering intensity from 
each magnetic domain; the data are shown for G = (— f>2>l) an d 3j f) m the K domain 
with Q = (5,-5,3)^ In this paper, we consider the K domain with Q = (3, 3, 3), not identical 
to the experimental one. In our K domain, the scattering vectors G = (—35 — 35 3) an d (3) 3> 3) 
correspond to the above experimental scattering vectors. Note that the (1,-1,0) crystal axis lies 
on the scattering plane at ip = in our definition of the geometrical factor (see the Appendix). 

Figure 4 (a) shows the calculated intensity as a function of photon energy for G = (— |, — 5, 3)- 
The azimuthal angle is set to be ip = —46°, such that the (1,0,0) crystal axis lies on the scattering 
plane. (This corresponds to the experimental geometry that the (0,1,0) crystal axis lies on the 
scattering plane.) For the a — > a' channel (lower panel), only one peak is found at the pre-iT-edge 
position, in agreement with the experiment. This peak arises from the quadrupole process, since 
the dipole process is forbidden. The amplitude from the quadrupole process interferes with the 
non-resonant term, leading to an anti-resonant dip in the spectra. For the a — > tt' channel (upper 
panel), a large peak is found at the main i^-edge position in addition to the pre-edge peak, in good 
agreement with the experimental curve after correction for absorption.!'© The intensity of the 
main peak depends on the S domains; for example, it is very small for domain S3. 

Figure 4 (b) shows the calculated intensity as a function of photon energy for G = (5, 3, §)■ The 
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Fig. 4. Magnetic scattering intensity as a function of photon energy (a) for ip = —46°, G = (— |, —h, |), and (b) 
for ip — 136°, G = (|, |, | ). The solid, dotted and broken lines represent the intensities from domains Si, S2, 
and 53 in the same K domain with Q = (i |, |). The zero of energy is at the main peak. The upper and lower 
panels show the intensity for the a — > ir' channel, and that for the a — > a' channel, respectively. The insets show 
the experimental curves of ref. 3. The curve in the inset on the top panel in (a) is after absorption correction. 



azimuthal angle is set to be ip = 136°, such that the (1, 0, 0) crystal axis lies on the scattering plane. 
(This correspond to the experimental situation that the (0, 1,0) crystal axis lies on the scattering 
vector.) For the a — » a' channel (lower panel), only the pre-edge peak is found, in agreement with 
the experiment. The anti-resonant dip is found below the pre-edge peak for domain S\. Such a dip 
is actually observed in the experiment (see the inset). For the a — > n' channel (upper panel), the 
main peak is found in addition to the pre-edge peak, but its intensity is much smaller than that 
for G = (— i, — f>f)- As already pointed out J!* the intensity of the main peak is the smallest for 
domain S\. Experimentally, the main peak is not clearly seen. Combining this fact to the spectral 
shape of the pre-edge peak, we speculate that the experimental spectrum comes from domain Si. 

Figure 5 shows the RXMS intensity for the a — > ir' channel in the energy region of main peak, 
which is calculated with neglecting the intraatomic exchange interaction between the Ap and 3d 
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Fig. 5. Magnetic scattering intensity for the a — > it' channel in the energy region of main peak, calculated with 
neglecting the intraatomic exchange interaction between the 4p and 3d states. The solid, dotted and broken lines 
represent the intensities from domains Si , S2 and S3 . 



states. The intensities are nearly half of those shown in Fig. 4, indicating that the contribution of 
the intraatomic exchange interaction is comparable to that of the p-d mixing to the 3d states of 
neighboring Co atoms. Note that the latter contribution increases according to roughly the square 
of the strength of the p-d mixing. It is highly probable that the effect of the intraatomic exchange 
interaction is suppressed by a band formation of the 4p states, since the energy splittings caused 
by the exchange interaction are averaged out at the A and B sublattices in the states with small 
momenta. For making clear this point, however, it is necessary to perform a band calculation, 
which has not benn completed. 

Figure 6 shows the calculated intensity as a function of the azimuthal angle for G = (— 5, — 5, |) 
and (5, 5, §)■ We find rather complicated behavior; for example, for the a — ► ir' channel at G = 
(5,5,!), the main peak grows considerably while the pre-edge peak shrinks, with increasing values 
of ip. One reason for such complicated behavior is that the scattering vector is not along the crystal 
axes. The experimental data are not available to make comparison. 

§5. Concluding Remarks 

We have analyzed the recent experiment for the RXMS intensity for CoO using a cluster model. 
The spins and orbitals are ordered with the same wavevector under strong spin-orbit coupling. 
The point different from the case of ordinary orbital order is that the scattering tensor becomes 
antisymmetric. The cluster calculation gives a reasonable description for the RXMS intensity. We 
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Fig. 6. Magnetic scattering intensity as a function of azimuthal angle ip; (a) for G = (— |,— 3,3), and (b) for 
G = (5, i, |). The upper panel is for the main K-edge (MKE) peak, and the middle and lower panels are for the 
pre if -edge (PKE) peak. The solid, dotted and broken lines represent the intensities from domains Si, S2 and S3. 



have found that the Ap states of Co are modified not only by the intraatomic exchange coupling to 
the 3d states but also by the orbital moment of the 3d states of neighboring Co atoms through the 
p-d mixing. Thereby the orbital moment is induced in the Ap states, and a considerable intensity is 
brought about through the dipole process. The p-d mixing mechanism reminds us of the magnetic 
circular dichroism (MCD) for the K-edge absorption in the ferromagnetic metals Fe, Co, and Ni,© 
where the MCD intensity comes mainly from the orbital moment induced in the Ap states through 
the mixing to the 3d states of neighboring Co atoms.!!© 

We have also found the pre-edge peak in the RXMS intensity. It originates from the quadrupole 
process. This contrasts with the pre-edge peak in the resonant scattering for the orbital order in 
LaMnOs, which is predicted to be generated by the dipole process.© We have also calculated the 
azimuthal angle dependence of the RXMS intensity, which shows rather complicated behavior. We 
are waiting for future experiments to make comparison. 

In this paper, we have used a simple cluster model for analyzing qualitatively the experimental 
data. Since the Ap states are considerably extended in space, it is much better to treat them as a 
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band. We speculate that the effect of the intraatomic exchange interaction is suppressed by a band 
formation of the 4p states. A study based on the band calculation is now under progress. 
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Appendix: Geometrical Factors 

We outline the derivation of the geometrical factors in the coordinate frame fixed to the crystal 
axes and summarize the result here. 

A.l Dipole Transition 

In the coordinate frame (x',y',z') that the z' axis is the propagation direction of photon, the 
dipole operators for circular polarization are given by 

1 /47r\ 1 / 2 

iy') = T (y) Y l±1 (9',<py, (A-l) 

where r', 9' , ip' are corresponding polar coordinates. We represent these operators in the coordinate 
frame (x" ,y" , z") that the z" axis is along the scattering vector G and y" axis is contained in the 
scattering plane (Fig. 3). This is done by the relation,!!) 

W(0V) =^2D^ m ,(a,P, 1 )Y em (e", ( p"), (A-2) 

m 

where I = 1, m! = ±1 with the Euler angles a, (3, 7 being a = vr/2, f3 = ir/2 + 9, and 7 = —tt/2 for 
incident photons, and being a = vr/2, (3 = tt/2 — 9, and 7 = —tt/2 for scattered photons. According 
to the rotation of the crystal by azimuthal angle ip around the scattering vector G, we replace 

Y lm (9", <p") — Y lm (9", tp")e im -*". (A-3) 

Next we then express Yi m (9" ,<p"ys in the coordinate frame (x'" ,y"' , z'") fixed to the crystal axes 
with the help of eq. ( |A-2| ), with a = Sir/A, j3 = arcsin >/2/51, 7 = 0, for G = (—5, —5, |), and 
with a = 7r/4, [3 = arcsin a/2/27, 7 = 0, for G = (|, |, |). Finally we change the expressions in the 
circularly polarized base into those in the linearly polarized base. Thus we obtain the geometrical 
factor for incident photon (/i = a or ir) as 

{P C7 ) 1 = ±(cos^cos^ + sinV)/v / 2, 
(P a ) 2 = ±(cos/3cos^ - smip)/V2, 
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(P a ) 3 = -sin/3cosV>, (A-4) 
(P 7T ) 1 = ± [sin6*(cos/3sin^ - cosV>) + cos 9 sin (3} /y/2, 
(P w ) 2 = ± [sin 6(cos (3 sin ip + cos t/0 + cos sin /3] / \/2, 

(P 71 ")^ = (- sin 6* sin /3 sin ^ + cos6>cos/3). (A-5) 

The upper and lower signs in eqs. (A-4) and (A-5) correspond to G = (— |, — |, |) and (|, |, §), 
respectively. We also obtain the geometrical factor for scattered photon P ,p! (// = a' or ir') from 
eqs. (A-4) and (A-5) with replacing sin9 by — sin 6. Note that the scattering plane for ip = 
contains the (1, —1,0) crystal axis in this definition of the coordinate frame. 

A. 2 Quadrupole Transition 

In the coordinate frame (x',y',z') that the z' axis is the propagation direction of photon, the 
quadrupole operators for circular polarization are given by 

11 /ir\ 1 / 2 

--j=z'{x' ± iy') = T [-) Y 2±1 (0>, iff). (A-6) 

We represent these operators in the coordinate frame fixed to the crystal axes by following the 
same procedure as for the dipole transition. We change the expressions in the circularly polarized 
base into those in the linearly polarized base. Thus we obtain the geometrical factor for incident 
photon (/J. = a or 7r) as 

(Q a )i = — cos cos (3 cos 2ip — sin 9 sin (3 sin ip, 

{Q a ) 2 = (y/S/2) cos 9 sin 2 (3 sin 2ip + ( V3/2) sin 9 sin 2(3 cosip, 

(Q a ) 3 = =p cos 9 sin /3(cos /3 sin 2^ + cos 2ip)/y/2 

=F sin (cos 2(3 cos V' — cos (3 sin VO / \/2, 
(Q CT ) 4 = =p cos 9 sin /3(cos /3 sin 2ip - cos 2^) /a/2 

=F sin 9 (cos 2/3 cos ip + cos /3 sin ip) / \[2, 
(Q CT ) 5 = cos6»[l - (l/2)sin 2 /3]sin2V> - (1/2) sin 9 sin 2(3 cos ip, (A-7) 
(Q 7r ) 1 = -(1/2) sin 29 cos /3 sin 2ip - cos 20 sin/? cos V>, 
(Q 7r ) 2 = - ( V3/4) sin 26* sin 2 /3 cos 2-0 - (a/3/2) cos 2(9 sin 2(3 sin V 

- (a/3/2) sin 20[1 - (3/2) sin 2 /3], 
(Q u ) 3 = ±(1/2a/2) sin 29 sin /3(cos /3 cos 2iP - sin 2^) 

± cos 20 (cos 2/3 sin ip + cos /3 cos ip) / \p2 

T (3/4A/2)sin26>sin2/3, 
(Q 7r ) 4 = ±(1/2 V2) sin 20 sin /3(cos /3 cos 2^ + sin 2ip) 

± cos 29 (cos 2/3 sin ^ — cos (3 cos ^) / \[2 
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T (3/4^2) sin 29 sin 2/3, 
(Q n ) 5 = -(l/2)sin20[l - (1/2) sin 2 0} cos 2^ 

+ (1/2) cos 29 sin 2/3 sin V - (3/4) sin 29 sin 2 /3, (A-8) 

where (3 = arcsin a/2/51 for G = (— i, — |, |), and /3 = arcsin a/2/27 for G = (|, |, |). The upper 
and lower signs in eqs. (A-7) and (A-8) correspond to G = (— |,— |, |) and (|, |, §). We also 
obtain the geometrical factor for scattered photon Q >tx ' (// = a' or n') from eqs. (A-7) and (A-8) 
with replacing sin6> and sin 29 by — sin 9 and — sin 29, respectively. 
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